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ABSTRACT 

We discuss a 175 deg 2 spectroscopic survey for blue horizontal branch (BHB) stars in the Galactic 
halo. We use the Two Micron All Sky Survey (2MASS) and the Sloan Digital Sky Survey (SDSS) 
to select BHB candidates, and find that the 2MASS and SDSS color-selection is 38% and 50% effi- 
cient, respectively, for BHB stars. Our samples include one likely run-away B7 star 6 kpc below the 
Galactic plane. The global properties of the BHB samples are consistent with membership in the 
halo population: the median metallicity is [Fe/H]=— 1.7, the velocity dispersion is 108 km s _1 , and 
the mean Galactic rotation of the BHB stars 3 < \z\ < 15 kpc is —4 ± 30 km s _1 . We discuss the 
theoretical basis of the Preston, Shectman & Beers My-color relation for BHB stars, and conclude 
that intrinsic shape of the BHB My-color relation results from the physics of stars on the horizontal 
branch. We calculate the luminosity function for the field BHB star samples using the Efstathiou, 
Ellis, & Peterson maximum- likelihood method which is unbiased by density variations. The field 
BHB luminosity function exhibits a steep rise at bright luminosities, a peak between 0.8 < My < 1.0, 
and a tail at faint luminosities. We compare the field BHB luminosity functions with the luminosity 
functions derived from sixteen different globular cluster BHBs. Kolmogorov-Smirnov tests suggest 
that field BHB stars and BHB stars in globular clusters share a common distribution of luminosities, 
with the exception of globular clusters with extended BHBs. 

Subject headings: stars: horizontal-branch — Galaxy: stellar content — Galaxy: halo 



1. INTRODUCTION 

Mapping the stellar halo requires objects that are suf- 
ficiently luminous to observe at large distan ces, yet com- 
mon e nough to sample the halo densely. In lBrown et al.l 
(2003), hereafter Paper I, we introduced the Century 
Survey Galactic Halo Project, a photometric and spec- 
troscopic survey from which we selected blue horizon- 
tal branch (BHB) stars as probes of the Milky Way 
halo. BHB stars meet our criteria for tracer sam- 
ples: they are intrinsically luminous and are quite nu- 
merous, with a number density in the halo that ex- 
ceeds that of RR Lvraes by roug hly a factor of ten 
ijPreston. Shectman. fc Beerall991lT . The spectral types 
of BHB stars are typically around AO, bluer than most 
competing stellar populations. As a result, candidate 
BHB stars in the halo are relatively easy to select by 
broadband colors alone. 

In Paper I we described the detailed stellar spectral 
analysis techniques developed for the Century Survey 
Galactic Halo Project. In this paper we investigate the 
mean Galactic rotation, metallicity, and luminosity func- 
tion of the halo BHB stars in the context of a comple- 



mentary 175 deg 2 spectroscopic survey. This new survey 
extends the work of the original Century Survey Galac- 
tic Halo Project by making use of two large-area, multi- 
passband imagi ng surveys: (1) th e Two Micron All Sky 
Survey r2MASS: ICu"tri et al.l2003[) and (2 ) the Sloan Dig- 
ital Sky Survey fSPSS: lYork "et^l27i00h . 

Previous spectroscopic surveys of field BHB stars 
have identified BHB stars over large ( several 10 3 
deg 2 ) areas of sky to shallower depths ijPierl Il983t 
iWilhelm etafl ll999b(L or over s mall (~10 2 deg 2 ) ar- 
eas of sky to greater depths JSOTnniCTjLarseiie^jl 
19891 lArnold fc Gilmorel 119921: [Kinman et alJ 1199 
Clewlev et al.ll2004t iKinman et al.ll2004[) thaiTthe Cen- 



tury Survey Galactic Halo Project. The exception is the 
recently published sample of 1170 BHB stars observed 
by the SDSS as mis-identifie d quasars or as fille r ob- 
jects in low density regions ijSirko et alJ l2004albf l. In 
comparison, our spectroscopic survey of BHB stars is 
cleanly selected and 100% complete within our color and 
magnitude selection limits. Combined with the original 
Century Survey sample, we have 157 spectroscopically- 
identified BHB stars over 239 deg 2 of sky. 
In §2 we describe the sample selection and spectro- 
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Fig. 1. — Sky map in Galactic coordinates. The Century Survey 
slice 8^5 < RA B i950 < 13 h 5, 29?0 < Dec B i950 < 30?0 is located 
in the north Galactic hemisphere. Our new survey slice 23 1 ?0 < 
RAj2000 < 3^67, -L°25° < Decj2ooo < +1?25 is located in the 
south Galactic hemisphere. 
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scopic observations of the new 175 deg 2 region and dis- 
cuss selection efficiencies for BHB stars. In §3 we discuss 
the basis of BHB luminosity-color-metallicity relations, 
and analyze the global kinematic and abundance prop- 
erties of our BHB samples. In §4 we calculate the lu- 
minosity functions for our field BHB star samples, and 
compare them with luminosity functions derived from 
globular cluster data. We summarize our results and 
conclude in §5. 

2. SAMPLE SELECTION 

2.1. Selection Region 

The original Century Survey Galactic Halo Project 
contains BHB stars selected by (V — R)o or (J — H) 
colors. Here we make use only of the sample selected 
with (V — R)n < 0.3, the "Century Survey" sample 
IjBrown et all2 0031. The original Century Survey sample 
covers a 1° x 64° slice located 8^5 < ai 950 < 13^5, 29° < 
^1950 < 30° and contains 39 spectroscopically-confirmcd 
BHB stars in the magnitude range 13 < Vq < 16.5 mag. 

Here we select BHB candidate stars from the 2MASS 
and SDSS surveys in a complementary region located 
along the celestial equator 23 h m s < a J2000 < 3 h 40 m s , 
-1°15'0" < <5j20oo < +1°15'0". Fig. □ is a plot of this 
70° x 2.5° region in Galactic coordinates. The survey 
is located predominantly at b < —45°, in a regio n that 
cleanly samples the halo in the lBrown et alJ (|2P04) BHB- 
candidate maps. 

2.2. 2MASS Selection 

The 2MASS catalog provide s uniform JHK pho tom- 
etry over the entire sky. In iBrown et alJ l)2004|) . we 
matched the original Century Survey sample to 2MASS 
and showed that 2MASS colors select A-type stars with 
~80% efficiency. The A-type stars are all good BHB 
candidates in our high Galactic latitude survey region. 

We have selected 90 BHB candidates from the 2MASS 
catalog in the magnitude range 12.5 < Jo < 15.5; BHB 
candidates have colors in the ranges —0.2 < (J — H)o < 
0.1 an d -0.1 < (H - K) < 0.1, following Brown et all 
(2004). Our upper color limits result in a high selec- 
tion efficiency but a reduced completeness for BHB stars. 
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Fig. 2.— The 2MASS-selected BHB candidate sample. Panel 
(a) shows the distribution of (J — H)o and Jo; the solid box shows 
the sample selection region. Panel (b) shows the distribution of 
(«' — g')o and (</ — r')rj colors; the dashed box shows the SDSS- 
sample selection region for comparison. 



Comparison with the original Century Survey sample 
shows that our color selection samples 65% of the BHB 
population IjBrown et al.H200"4j) . 

It is important to note that we have selected objects us- 
ing d e-reddened colors and magnitudes, us ing extinctions 
from lSchleeel. Finkbeiner. fc Davisl l)1998j) . The surface 
density of the 2MASS-selected BHB candidates is 0.5 
deg -2 . 

We have matched up our 2MASS-selected BHB can- 
didates with the publicly available SDSS data: SDSS 
photometry presently exists for 65 of the 90 objects. Ap- 
proximately half of the matched objects have SDSS colors 
consistent with early A-type stars; the remainder follow 
the stellar locus to F-type stars (see Fig- El ■ 
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2.3. SDSS Selection 

The SDSS has released five passband photometry for 
limited areas of the sky that can be used to select A-type 
stars efficiently. We selected 194 BHB candidates in the 
magnitude ra nge 15 < q' n < 17 fro m the SDSS Early 
D ata Release Ipt oughto n et al.ll2002|) and Data Release 
1 (jAbazaiian et al.ll2003j) . We follow lYannv et all (|2000f) 
and select BHB candidates with —0.3 < (g' — r')o < 0.0 
and 0.8 < (u' - g') < 1.5. BHB candidates that fall 
outside the selection box in Fig. were objects originally 
selected by "model" magnitudes from the Early Data 
Release; here we plot Data Release 1 Petrosian magni- 
tudes that we find are better behaved at bright magni- 
tudes. The surface density of the SDSS-selected BHB 
candidates is 1 deg~ 2 . There is no overlap of these ob- 
jects with the 2MASS-selected sample, even though both 
samples cover the same region of sky. 

We looked up available 2MASS photometry for the 
SDSS-selected BHB candidates and found matches for 
188 of the 194 objects. A handful of SDSS stars sat- 
isfy the 2MASS-selection in Jo and (J — H)q but are 
rejected by (H — K)q. Thus the lack of overlap between 
the 2MASS- and SDSS-selected samples is likely due to 
the extreme uncertainties in 2MASS colors for the fainter 
16th and 17th magnitude SDSS stars (see Fig. • Inter- 
estingly, errant G-type stars found in the SDSS sample 
are cleanly identified by 2MASS photometry as bright 
and red (J — H)q ~ 0.35 stars. This comparison sug- 
gests that some bright (15th to 16th magnitude) SDSS 
stars are likely saturated, and thus have erroneous re- 
ported magnitudes. The on-line documentation for the 
SDSS data archive now describes a series of flags that 
can be used to avoid such saturated objects. 

To understand our completeness for BHB stars requires 
a better understanding of the SDSS saturation prob- 
lem. We start by selecting all stars with A-type col- 
ors along the celestial equator in SDSS Data Release 2 
ijAbazaiian et al.ll2004|) . We find that saturated objects 
have discrepant fr' — i')o colors for A-type stars. The 
solid line in Fig.@]shows the fraction of objects with dis- 
crepant (r' — i')o > 0.3. We then re-select all A-colored 
stars, but this time using the photometry flags to select 
objects only with clean photometry. The dashed line in 
Fig- El shows the fraction of objects with clean photome- 
try. 

Figure 0] shows that half of all A-colored stars with 
15 < g' Q < 15.5 have erroneous photometry and are not 
A stars at all. For these objects to have A-type colors 
in (u' — g') and (g' — r') but not in (r' — i') suggests 
that the g' band is saturated. This result also suggests 
that half of the real A-colored stars may be missing in 
this magnitude range. Selecting for clean photometry 
removes the erroneous objects, but may also reduce the 
completeness of the sample. The fraction of discrepant 
A-colored stars drops to ~10% at g' = 15.75 (see Fig.0J, 
and the clean photometry selection maintains this level 
of apparent incompleteness to g' Q = 17. 

2.4. Spectroscopic Observations 

During the fall 2003 observing season we obtained a 
spectrum for each BHB candidate in the 2MASS- and 
SDSS-selected samples. Spectroscopic o bservations were 
obtained with the FAST spectrograph (Fabr icant et alJ 
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Fig. 3.— Same as Figure EJ but for the SDSS-selected BHB 
candidate sample. Objects outside the SDSS-sample selection re- 
gion (solid box) are BHB candidates originally selected by "model" 
magnitudes from the Early Data Release; here we plot Data Re- 
lease 1 Petrosian magnitudes that we find are better behaved at 
bright magnitudes. 



1998) on the Whipple 1.5 m Tillinghast telescope. We 
used a 600 line mm -1 grating and a 2 arcsec slit to obtain 
a resolution of 2.3 A and a spectral coverage from 3400 to 
5400 A. Typical signal-to-noise (S/N) ratios were 30/1 
in the continuum for objects brighter than 16th magni- 
tude, decreasing to S/N—l 5/1 for the 17th magnitude 
objects. This S/N is adequate to measure the Balmcr 
lines and the Balmer jump - the primary surface-gravity 
indicators we employ for BHB stars. Paper I contains 
details of the data reduction. We measure spectral types 
and radial velocities, and derive metallicities, effective 
temperatures, and surface gravities, from the spectra of 
the total sample of 284 objects. 
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Fig. 4. — Fraction of all A-colored stars in the SDSS DR2 equa- 
torial region with bad photometry flags (dashed line) or discrepant 
(r' — i')o > 0.3 color (solid line). 
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SELECTION EFFICIENCIES 
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Fig. 5. — The four BHB classification techniqu es applied to 
our sample: (a) the modifie d ll^nman et alj tl99l) method, (b) 
thelWilhelm. Beers. &; Gray ( 1999a) metho d, (c) the l C lewlev et alj 
120021) On i^-Color method, and (d) the lClewlev et alj I2002I) Scale 
width-Shape method. Solid circles mark the BHB stars; open cir- 
cles mark the high surface-gravity A-type stars. 



2.5. BHB Classification 

The major difficulty in using BHB stars as probes 
of Galactic structure is the need to distinguish re- 
liably between low surface-gravity BHB stars and 
higher surface-gravity A-type dwarfs and blue strag- 
glers. Although investigators once thought blue 
stragglers were a minor c omponent of the halo 
population, recent studies l | No rris fc Hawkins] |1991t 
Preston. Beers, fc Shectmadll994IWilhelm et alll 999b: 
Brown et alJ l2003t iClewlev et alJ l2004[) demonstrate 
that a surprisingly large fraction of faint stars in 
the color range associated with BHB stars are indeed 
high-gravity stars, many of which are blu e stragglers 
ijPreston fc Snedenll200alCarnev et alJl2005|) . 

Our classification of BHB stars is identical to the ap- 
proach d escribed in Paper I. I n brief, we apply the tech- 
niques of lKinman et alj l|1994l) , iWilhelm. Beers, fc Gravl 
jl999a|) . and IClewlev et all l)2002j) to identify low 
surface-gravity BHB stars. We identify objects that sat- 
isfy three or more of the four classification techniques as 
BHB stars (see Fig.|5J). We find a total of 118 BHB stars 
across our 175 deg 2 survey region. 

2.6. Sample Selection Efficiencies 

Table ^ summarizes sample selection efficiencies. The 
2MASS-selected sample contains 34 BHB stars (out of 
90 candidates) for a net selection efficiency of 38%. The 
total number of A-type stars is about twice the number 
of BHB stars, or 78% of the 2MASS-selected sample. Of 
the remaining non A-type objects, 7% of the stars in 
the 2MASS-selected sample are B-type stars; 15% of the 
stars in this sample are F-type stars. 

The SDSS-selected sample contains 84 BHB stars (out 



of 167 candidates) for a net selection efficiency of 50%. 
We ignore the 27 G-type stars in this calculation as these 
stars can presumably be rejected by saturation flags. The 
total number of A-type stars is about twice the number of 
BHB stars, or 92% of the SDSS-selected sample (exclud- 
ing the G-types). Of the remaining non A-type objects, 
4% of the stars in the SDSS-selected sample are B-type 
st ars; 4% of the stars in this sample are F-type stars. 

Sirk o et aT] l|2004aj) have recently published a "strin- 
gent" color selection for BHB stars. Applying the strin- 
gent color cut to our full SDSS-selected sample would 
yield 55 BHB stars selected from 81 candidates for a net 
selection efficiency of 68%, but a completeness of only 
65% compared to the full SDSS-selected sample. 

2.7. Unusual Objects 

In Paper I we identified a number of unusual objects, 
including white dwarfs, subdwarfs, and B-type stars, 
within our survey of blue stars in the halo. The 2MASS- 
and SDSS-selected samples, by comparison, contain a 
handful of B stars, but do not include any white dwarfs or 
subdwarfs. The lack of white dwarfs may be explained 
by the more restrictive color selection we used for the 
2MASS- and SDSS-selected samples. Moreover, the B- 
type stars in the 2MASS- and SDSS-selected samples are 
almost entirely late B8 and B9 stars. These late B-types 
are potentially all hot horizontal-branch stars, but are 
very difficult to classify because the horizontal branch 
crosses the main sequence at this location in the H-R 
diagram. 

The earliest B-type star in our samples is CHSS 1645, 
classified as B7. As the earliest B-type star in our sam- 
ples, CHSS 1645 is the most likely object to be a true 
B star rather than a hot horizontal-branch star. Assum- 
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ingCHSS 1645 has solar metallic ity, with My ~ —0.6 
llCoxll2000l) and (V - J) = -0.3 l|Kenvon fc Hartmannl 
1995), we estimate that it is located 6 kpc below the 
Galactic plane. This places CHSS 1645 among the class 
of stars known as "run-away B-type" stars. The star 
CHSS 1645 is located at b = -60°, hence its +73 km 
s _1 radial velocity points predominantly in the negative 
z direction perpendicular to the plane of the Galaxy. If 
its radial velocity is the majority of its full space motion, 
it takes 10 8 years for CHSS 1645 to trave l 6 kpc from the 
Galactic plane. A B7 star has ~4 M^ llCoxl 120001) and 
a lifetime ~2xl0 8 yr (.Bo wers fc Deeminelll984l) Thus 
CHSS 1645, a likely run-away B7 star, has a lifetime con- 
sistent with its travel time from the disk. 

3. GLOBAL PROPERTIES 

To map the Galactic halo requires knowing the intrinsic 
luminosities of BHB stars. BHB stars are standard can- 
dles with luminosities that depend on effective tempera- 
ture (color) as well as metallicity. We begin by discussing 
the physical basis of the BHB luminosity dependence on 
color (Section 3.1). We then present the observed distri- 
bution of metallicities derived from our spectra (Section 
3.2). Using our colors and metallicities, we compute in- 
trinsic luminosities for our field BHB stars and investi- 
gate their spatial distribution (Section 3.3). Finally, we 
investigate the mean Galactic rotation of our halo sam- 
ples (Section 3.4). 

3.1. BHB Luminosity- Color Dependence 

BHB stars share a common physical origin. They are 
stars that have evolved off the red giant branch and are 
burning helium in their cores with a hydrogen burning 
shell. The bolometric luminosity of a BHB star depends 
on t he core mass, the stell ar mass, and the metallicity 
(e.g. iDemaraue et al.l 120001. More massive BHB stars 
have larger hydrogen-rich envelopes and are cooler than 
less massive BHB stars. The variation of effective tem- 
perature with stellar mass yields a robust relation be- 
tween optical luminosity and B-V color: bluer BHB stars 
a re fainter than red BHB stars . 

iPreston. Shectman. fc Beers! (|1991|) . hereafter PSB, 
provide an empirical BHB luminosity-color relation from 
a fit to fifteen globular cluster BHBs. Fig. shows the 
PSB My-color relation as a solid line. The relation is 
normalized to [Fe/H]=— 2.3 using the My-metallicity re- 
lation described below. Globular clusters exhibit a wide 
range of BHB morphologies, evident in the 0.25 mag scat- 
ter of the points in PSB's Fig. 5. Because the large scat- 
ter may result from the physics of globular clusters rather 
than from the physics of BHB stars, we next consider a 
theoretical My-color relation. The theoretical My-color 
relation provides a physical basis for the empirical PSB 
relation. 

To construct a theoretical Afy-color relation, we 
use the horizon t al br anch evolutionary tracks from 
IDemaraue et all (J2000) . For models with Z = 10~ 4 
(equivalent to [Fe/H]= —2.3) we adopt the luminosity 
at time zero and derive colo rs and bolometric tempera- 
tures from published tables llKen von fc Hartmannl 119951 
iGreen fc Demarauelll996l iLeieune et al .11 19981) . The re- 
sulting My-color relation for Z = 10 -4 is the dashed 
line in Fig. The theoretical My-color relation is re- 
markably similar in shape to the empirical PSB relation 
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from globular clusters , the dashed line shows the theoretical rela- 
tion derived from the Demarque et all 120001) stellar models, and 
the dotted line shows our toy model of bolometric corrections 
(Kcnvon & Hartmann 1995) for a constant luminosity star. 



in the (B — V)q > —0.1 region covered by our BHB star 
samples. 

As a consistency check, we plot a third line in Fig. 
El that is simply the bolometric correction for a star 
with constant luminosity. We use the bolometric correc- 
tions f or main sequence stars from lKenvon fc Hartmanrl 
(1995), and add 0.5 mag to match the bolometric cor- 
rection to the empirical and theoretical My-color re- 
lati ons at the red end. Inter estingly, the shape of 
the iKenvon fc Hartmannl (|l995) bolometric corrections 
is similar to both the empirical and theoretical BHB 
My-color relations, except that the slope of the bolo- 
metric correction curve is too shallow at the blue end. 
We expect this systematic difference because blue BHB 
stars are intrinsically less luminous than red BHB stars; 
we have assumed a co nstant luminosity s t ar. Bolo- 
metric correct i ons fr om IGreen fc Demaraud (|1996j) and 
ILeieune et alJ l|1998D yield similar results, with a typical 
scatter of 0.1 - 0.2 mag. This toy model shows that the 
primary ingredient in the BHB My-color relation is the 
bolometric correction for BHB stars. 

Thus the physics common to all BHB stars leads to 
a general BHB My-color relation, albeit with an intrin- 
sic spread resulting from age and metallicity. The My- 
color relation depends on age because the luminosity and 
effective temperature of a BHB star evolve with time. 
The My-color relation has a wel l know n dependence on 
metallicity, but IDemaraue et al.l l)2000t) argue for an ad- 
ditional spread in My at a given metallicity due to BHB 
morphology. The morphology effect is strongest for a 
metal poor [Fe/H]< —2 BHB with blue morpho l ogy (H B 
type index=+l). According to Dcma roue et alJ (|2()00), a 
metal poor, blue BHB is actually ~0.1 mag brighter than 
the standard luminosity-metallicity relation predicts. We 
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Fig. 7. — Distribution of BHB [Fo/H] for (a) the Century Survey 
sample, (b) the 2MASS-selected sample, and (c) the SDSS-selected 
sample. Dashed lines indicate the median [Fe/H], The error in 
[Fe/H] is 0.25 dex, the same size as our bins. 



conclude that the BHB My-color relation has an intrin- 
sic shape due to the physics of the horizontal branch, 
with an intrinsic spread of 0.1 - 0.2 mag. For purposes of 
discussion, we use the empirical PSB My-color relation 
to estimate BHB luminosities in the rest of this section. 

3.2. Metallicities 

We measure metallicities for BHB stars as described 
in Paper I. We use three diffe rent techniques: the 
line indices of iBeers et alJ l)1999j) . the equivalent width 
of Ca II K plus a chi-square comparison between 
metallic-line regions in syntheti c and observed spectra 
(Wilhclm. Bee rs, fc Oravl Il99 9a ? ) . and a n optimization 
meth od that fits the entire spectrum ijAllende Prietol 
2002§. The three techniques are in good a greement 
with 0.25 dex uncertainty l|Brown et all 120031) . The fi- 
nal metallicity is the average of the three techniques; we 
adopt ±0.25 dex as the error in the final metallicity. 

Figure[7|plots the observed distribution of metallicities 
we measure in the original Century Survey sample, the 
2MASS-selected sample, and the SDSS-selected sample. 
The median [Fe/H] of the BHB samples are indicated by 
the dashed lines in Fig. and range from [Fe/H] =—1.47 
to —1.75. Our metallicity- measuring techniques are lim- 
ited to the range —3 <[Fe/H]< 0, and so peaks in Fig. 
0at [Fe/H]=— 3 and are likely stars with lower/higher 
metallicities piling up at the limits. We caution that the 
distributions in Fig. [7] are the observed and not the in- 
trinsic distributions: because BHB luminosities have a 
mild dependence on metallicity, stars of different [Fe/H] 
are sampled from different volumes of space (see below). 
However, it is clear that our BHB samples are predom- 
inantly metal-poor and therefore consistent with a halo 
population. 

A Kolmogorov-Smirnov (K-S) two-sample test pro- 
vides a simple way to evaluate whether the different BHB 
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Fig. 8.— Distribution of [Fe/H] versus z for (a) the 2MASS- 
and (b) the SDSS-selected samples. The dotted lines represent the 
magnitude limits for a BHB star at the median Galactic latitude 
b = —57° and at the median color BV0 = 0.03 of the samples. 
Some stars fall beyond the "limits" because there is a spread of 
Galactic latitude and color in our samples. 



samples are drawn from a common parent distribution 
of [Fe/H]. The K-S two-sample test works by sorting an 
observed quantity, such as metallicity, and then compar- 
ing the cumulative distributions of two different samples 
with one another. The likelihood is calculated for the null 
hypothesis that the two distributions are drawn from the 
same parent distribution. We test the metallicity distri- 
butions in Fig. din a pairwise fashion, and find likelihood 
values ranging from 14% to 60%. Thus our BHB samples 
are consistent with the null hypothesis that the metallic- 
ity distributions come from the same (halo) population. 

3.3. Spatial Distribution 

We calculate luminosities for our fie ld BHB stars 
using the My {BHB) relation from iClewlev et alJ 
(2004). This relation assumes the Hipparcos-devwed 
zero p oint, My(RR) = 0.7 7 ± 0.13 at [Fe/H] = 
-1.60 (jGould fc Popowskil 19981) . a My-metallicity slope 

0. 214±0.047 based on RR Lyrae star s in the Large Magel- 
lanic Cloud l|Clementini et al.ll200^ . and the PSB cubic 
relation in (B — V)o to provide the temperature correc- 
tion. Although the PSB Afy-color relation was derived 
for globular cluster BHB stars, the shape of the rela- 
tion reflects the physics common to all BHB stars, as 
explained above. Note that we do not measure (B — V)q 
directly. For the SDSS sample, we are able to make ac- 
curate estimates of (B — V)q from SDSS colors. For the 
2MASS sample, we use 2MASS photometry and Balmer 
line strengths to estimate (B — V)o as described in Paper 

1. We refer to these (B — V)q estimates as BV0. From the 
derived luminosities we compute distances. We expect 
the relative distances of our BHB stars have a precision 
of -6%. 

Because the luminosity of a BHB star is dependent on 
metallicity, the depths reached by our flux-limited sam- 
ples are dependent on metallicity. Fig. |H1 shows the dis- 
tribution of [Fe/H] for the 2MASS- and SDSS-selected 
samples as a function of z, the distance above or be- 
low the Galactic plane. The dotted lines represent the 
magnitude limits for a BHB star at the median Galactic 
latitude b = —57° and at the median color BV0 = 0.03 
of our samples. Metal-rich BHB stars are intrinsically 
fainter than metal-poor BHB stars. Fig.[H]shows that we 
sample BHB stars with [Fe/H]= -1 to 82% the depth of 
BHB stars with [Fe/H]= -3. 
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Fig. 9.— Distribution of BVO versus z for (a) the 2MASS- 
and (b) the SDSS-selected samples. The dotted lines represent the 
magnitude limits for a BHB star at the median Galactic latitude 
b = —57° and at the median metallicity [Fe/H]= — 1.7 of the sam- 
ples. Some stars fall beyond the "limits" because there is a spread 
of Galactic latitude and metallicity in our samples. 



The BHB luminosity dependence on color is stronger 
than the dependence on metallicity. Thus there is a 
strong selection bias with stellar color. In Fig. we plot 
the distribution of BVO color for the 2MASS- and SDSS- 
selected samples as a function of z. The dotted lines rep- 
resent the magnitude limits for a BHB star at the median 
Galactic latitude b = —57° and at the median metallic- 
ity [Fe/H]= —1.7 of our samples. Fig. [3] shows that we 
sample BHB stars with BVO = -0.1 to only 64% of the 
depth that we detect BHB stars with BVO = +0.1. The 
intrinsically faintest BHB stars are the bluest BHB stars 
hooking down off the horizontal branch in an H-R dia- 
gram. These faint BHB stars are sampled in a smaller 
volume than the more luminous BHB stars in our sam- 
ples. 

We now plot the spatial distribution of the original 
Century Survey the 2MASS-selected, and the SDSS- 
selected BHB samples (Fig. I10|l. Spatial distribution is 
traditionally displayed in a wedge plot for survey slices 
like ours. However, a wedge plot is inappropriate in the 
context of the Galaxy where a slice in celestial coordi- 
nates cuts across varying Galactic latitudes. The density 
of halo and disk populations is a strong function of both 
R, the distance along the Galactic plane, and z (e.g. 
iSiegel et all 120021) . Thus in Fig. EH we plot the distri- 
bution of BHB stars as a function of R and z. Indeed, 
the observed distribution of BHB stars clearly depends 
on both R and z. 

The 2MASS- and SDSS-selected BHB candidates sam- 
ple complementary ranges of distances, but the overlap is 
unfortunately minimal. The 2MASS and SDSS catalogs 
formally overlap between 15 and 15.5 magnitude. How- 
ever, in this magnitude range, the 2MASS BHB selection 
suffers from incompleteness due to large color errors and 
the SDSS BHB selection likely suffers from incomplete- 
ness due to saturation problems. 

3.4. Mean Galactic Rotation 

There is wide disagreement in the literature on 
whether the stellar halo rotates significantly. Pre- 
vious surveys have found evidence for (1) no 
halo r otation dLaaxlen et aU Il99fjt iGou ld fc Po powskil 
19981 iMartin fc Morrison! 119981: iGilmore et all 120021: 
" Sirko et alJ l2004bD. (2) a small prograde rotation 
llChiba fc Beersl I2000D . and (3) retr ograde rotation 
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Fig. 10. — Distribution of BHB stars in distance from the Galac- 
tic center along the Galactic plane, R, and distance above/below 
the Galactic plane, z. Panel (a) shows the BHB stars in the orig- 
inal Century Survey slice. Panel (b) shows the BHB stars in the 
2MASS-selected (triangles) and SDSS-selected (crosses) samples. 
The solid lines indicate heliocentric distance limits of 2 kpc and 
(a) 14 kpc or (b) 17 kpc. 



(IMaiewskil 119921: IMaiewski et all 119961 ISpaena et all 
120031 IKinman et alJl2004j) . Interestingly, all the mea- 
surements of retrograde rotation come from surveys of 
the north Galactic pole. By comparison, the measure- 
ments of no rotation come from surveys covering many 
directions in the sky. 

Our BHB samples cover a wide range of Galactic lati- 
tude and longitude, and so provide us with a "fair" sam- 
ple of the halo. Although high Galactic latitude stars are 
not ideal for measuring the rotation of the stellar halo, 
the 2MASS- and SDSS-selected samples include a num- 
ber of stars near I ~ 90° that are sensitive to a systematic 
rotation of the halo. 

Figure ^2 shows the mean rotation velocity of the 
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Fig. 11— Galactic rotation velocity of the 2MASS- and SDSS- 
selected BHB stars, binned by distance below the Galactic plane. 



2MASS- and SDSS-selected BHB stars as a function of 
z. Each bin includes ~25 BHB stars that are first cor- 
recte d to the local standard of rest ijDehnen fc Binnevl 
19981) and then eva luated by the methodology of 
Frenk fc White! (|198flfl . We assume a solar rotation of 
220 km s _1 . We expect some contamination from the 
thick disk in the nearest BHB stars, and indeed there 
is a hint of prograde rotation in the 1 < \z\ < 3 kpc 
bin. The mean rotation velocities of the 3 < \z\ < 15 
kpc stars are, however, consistent with no rotation. The 
rotation velocity of the combined sample of 2MASS- and 
SDSS-selected BHB stars is -4 ± 30 km s" 1 . The ve- 
locity dispersion of the BHB stars is 108 km s , also 
consistent with a halo population. 

The BHB stars from the original Century Survey sam- 
ple cover a similar range of z. However, the Century 
Survey BHB stars are located towards the Galactic anti- 
center I ~ 200° and towards the north Galactic pole 
b > 60°. Thus the Century Survey stars provide very 
little leverage on halo rotation. When we include the 
Century Survey stars in the mean rotation velocity cal- 
culation, we find they add ~ 10 km s _1 of retrograde 
rotation to the bins. We conclude the mean rotation 
velocities remain fully consistent with no halo rotation 
within their errors. 



4. THE LUMINOSITY FUNCTION OF BHB STARS 

Knowledge of the intrinsic distribution of luminosities 
of field BHB stars is important for interpreting maps of 
the Galactic halo. Knowledge of the luminosity function 
is also important for understanding the intrinsic proper- 
ties of field BHB stars that cover a broad range of ob- 
served magnitude, color, and metallicity. The luminosity 
function describes the number of stars per unit volume 
in the luminosity interval L to L + dL. We describe the 
method we use to calculate the BHB luminosity function 
(Section 4.1), and discuss the role of the My-color rela- 
tion in our result (Section 4.2). We compare the luminos- 
ity function we determine for our field BHB stars (Sec- 
tion 4.3) with luminosity functions derived from globular 
clusters with BHBs (Section 4.4). 

4.1. Calculating the Luminosity Function 

We calculate the luminosity function of our field BHB 
stars using the non-parametric maximum-likelihood 
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Fig. 12. — We recover the intrinsic luminosity functions (lower 
panel) for four simulated intrinsic color distributions (upper panel), 
in this case using the PSB M^-color relation. The distribution of 
colors is the crucial element for the luminosity function. 



metho d of lEfstathiou et~aT] (|1988j) . The lEfstathiou et alJ 
( 1988) maximum-likelihood method is commonly used to 
calculate the luminosity function of galaxies in galaxy 
redshift surveys. We now apply this method to our sur- 
vey of BHB stars in the Galactic halo. The method does 
not simply count the numbers of stars at different lumi- 
nosities, but weights the contribution of each star by the 
relative volume in which it can be observed in a flux- 
limited sample. Specifically, the probability of a star at 
distance d falling into the luminosity range [L,L + dL] is 
equal to the luminosity function at L divided by the num- 
ber density of stars one expects to see in a flux-limited 
survey at distance d. The maximum-likelihood method 
works by maximizing the sum of these probabilities and 
solving for the best-fitting luminosity function. 

The density terms drop out in the maximum-likelihood 
formalism with two notable consequences. First, the 
maximum-likelihood method is unbiased by systematic 
density variations. The maximum-likelihood method 
does not require knowledge of the halo density distribu- 
tion p(R, z); it only requires that the luminosity function 
is independent of position in the sampled volume. Sec- 
ondly, the absolute normalization of the luminosity func- 
tion is lost and requires a separate computation. Because 
stellar density varies with position in the Milky Way and 
because our samples are too sparse to fit the density pro- 
file directly, we compute only the form of the luminosity 
function, and arbitrarily normalize the luminosity func- 
tions to unity. 

4.2. The Role of the My-color Relation 

The My-color relations (Fig. [fJJ) specify only how a 
particular color maps to a particular My. The distri- 
bution of colors is not at all specified by the My-color 
relation. The distribution of colors is the crucial element 
of the luminosity function. To illustrate this point, Fig. 
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Fig. 13. — Luminosity functions of halo BHB stars in (a) the 
Century Survey sample, (b) the 2MASS-selected sample, and (c) 
the SDSS-selected sample. The normalization is scaled so that the 
areas under the curves are equal to one. 



^| plots luminosity functions calculated for the following 
four intrinsic color distributions: (1) an uniform color 
distribution, (2) a Gaussian color distribution centered 
at (B — V)o=0.05 with a — 0.05 mag, and (3) exponential 
color distributions with scale length 0.1 mag peaking in 
the red and (4) the blue. Each model color distribution 
contains 1000 objects. For purposes of this calculation, 
we derive the intrinsic BHB luminosities using the PSB 
My-color relation, though our results are nearly identical 
for the other My-color relations in Fig. The bottom 
panel of Fig. 12 shows the luminosity functions resulting 
from the four color distributions. 

It is clear from Fig. ^] that the BHB luminosity func- 
tion depends dramatically on the distribution of BHB 
colors. Each simulated luminosity function in Fig.ll2lh.as 
a different shape, some with narrow distributions, others 
with long tails extending to faint luminosities. Moreover, 
the characteristic peaks of the luminosity functions vary 
in luminosity and total number of stars. 

Even though colors are the primary indicator of BHB 
luminosity, we cannot compare raw distributions of col- 
ors because our field BHB stars have different luminosi- 
ties and thus sample different volumes of space. To de- 
rive intrinsic properties requires knowing the luminosity 
function of our field BHB stars. 

4.3. The Field BHB Luminosity Function 

Figure 1131 shows the luminosity function of the origi- 
nal Century Survey, the 2MASS-selected, and the SDSS- 
selected samples, determined from the observations of 
color, metallicity, and apparent magnitude. We use 0.2 
magnitude wide bins to encompass any uncertainties in 
the My derivati on. Luminosit i es are derived with the 
previously stated IClewlev et al.l 1)20041) relation that uses 
the PSB My-color relation. All three luminosity func- 
tions display the same general shape: a steep rise at 
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Fig. 14. — Average luminosity function of our three halo BHB 
star samples, calculated using (1) the empi rical PSB My-co lor re- 
lation (solid line), (2) the theoretical IDemarque et alj 12000) BHB 
model s (dashed line), and (3) our toy model of bolometric correc- 
tions iKcnvon & Hartmann 1995) for a constant luminosity BHB 
star. 
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Fig. 15. — Apparent magnitude distribution of BHB stars in (a) 
the Century Survey sample, (b) the 2MASS-selected sample, and 
(c) the SDSS-sclcctcd sample. Extinction-corrected magnitudes 
are (a) Vq, (b) Jo, and (c) g' . Solid lines indicate the number of 
BHB stars predicted by our derived luminosity functions assuming 
a r -2 ' 5 halo density profile, normalized to the number of stars in 
each sample. 



bright luminosities, a peak between 0.8 < My < 1.0, 
and a tail at faint luminosities. 

We perform K-S tests to measure the likelihood that 
our field BHB luminosity functions are drawn from the 
same parent population. The K-S test applies to un- 
binned distributions; we thus multiply the luminosity 
functions (Fig. 113(1 by the number of objects in the sam- 
ples and distribute the My's uniformly across each bin. 
The resulting likelihoods range from 37% to 57%, sug- 
gesting that our BHB samples share a common parent 
population. These likelihoods also mean that the bi- 
modal distribution of luminosities in the Century Survey 
sample fFig.H3fa) is not statistically significant. Interest- 



10 



W. R. Brown et al. 



ingly, all three samples share the same median My ~ 1.0 
mag. This agreement is rather remarkable, given the in- 
dependent photometry of the three samples, and suggests 
that the Century Survey, 2MASS, and SDSS have con- 
sistent photometry. 

Because the K-S tests suggest our three BHB samples 
are drawn from the same halo population, we average 
the three BHB samples to obtain a more robust measure 
of the field BHB luminosity function. We multiply each 
luminosity function in Fig. 1131 by the total number of 
BHB stars in each sample, sum the luminosity functions, 
and then divide the result by the grand total of BHB 
stars. The result is plotted as the solid histogram in Fig. 

We re-compute the luminosity functions using the My- 
color re lations derived from ( 1 ) IDemaraue et all (2000) 
and (2) iKenvon k, Hartmannl l)1995|) . We show the re- 
sults as (1) dashed and (2) dotted histograms in Fig. 
1141 Because we are interested in the shape of the 
luminosity function, we adjust the zero-points of the 
My- color relations derived from | Demaraue et al.l (2000) 
and IKenvon fc Hartmannl l)1995f) to match the zero- 
point of the PSB relat ion. Interestingly, the theoretical 
IDemaraue et all l|2000|) and empirical PSB curves have 
very similar shape: a K-S test gives a 99% likelihood for 
the two samples to share a common distribution. Thus 
the shapes of the theoretical and empirical My-color re- 
lations are similar enough to have no apparent effect on 
the shape of the final BHB luminosity function. 

Knowledge of the BHB luminosity function allows us, 
in theory, to solve for the BHB density distribution. In 
practice, our relatively sparse samples do not provide an 
adequate constraint (see Fig. lT5|l . We note t hat the halo 
power laws and scale lengths published in iSiegel et al.1 
(2002) yield reduced \ 2 ~ 1 an d so appear consistent 
with the distribution of our field BHB stars. 

To check of the veracity of our average BHB luminosity 
function, we use the luminosity function in Fig. 1141 (solid 
histogram) to calculate the expected apparent magni- 
tude distributions of our BHB samples. Fig. El plots 
the observed number of BHB stars in the original Cen- 
tury Survey, the 2MASS-selected, and the SDSS-selected 
samples versus extinction-corrected apparent magnitude. 
The errorbars indicate y/~N uncertainties. The solid lines 
in Fig. El are the number of BHB stars predicted by the 
BHB luminosit y function, assum ing a r~ 2 ' 5 power-law 
density profile ijSiegel et al J 120021) . To set the normal- 
ization, we scale the predictions to the observed number 
of stars in each sample. There is good agreement in the 
predicted shape of the magnitude distribution and the 
observations. The one exception is the final 2MASS bin 
with 15 < Jo < 15.5. We attribute the observed under- 
density to larger photometric errors at faint magnitudes 
(see Fig. 121 that scatter BHB stars out of the narrow 
2MASS color-selection box, thereby reducing our com- 
pleteness. 

4.4. Comparison with Globular Clusters 

An additional insight into our field BHB luminosity 
function is provided by comparison with globular clus- 
ter data. The purpose of this comparison is not to sug- 
gest that the halo is made of disrupted globular clusters. 
Rather, because all BHB stars share a common physical 
basis, we inquire whether they exhibit a common par- 



ent distribution of luminosities. Globular cluster BHB 
morphologies are known to vary widely because of dif- 
ferences in metallicity, main sequence turn-off mass, and 
"second-parameter" effects. We expect that our wide 
area surveys of the halo will sample BHB stars from the 
full range of BHB morphologies. We now test whether 
field and globular cluster BHB stars share a similar or 
different distribution of luminosities by comparing the 
shapes and median My's of the BHB luminosity func- 
tions. 

In a brief example of the virtual observatory in ac- 
tion, we used the NASA ADS system ijKurtz et al.ll2000j) 
to locate suitable globular cluster data for comparison 
with our field BHB samples. We used the query "globu- 
lar cluster color magnitude diagram" and required that 
there be on-line data associated with the paper. The 
first (most rece nt) paper wh i ch me t all our requirements 
is the study bv lHargis et all ll200l of M12 (NGC 6218). 
We followed the data lin k to the CDS/ VIZIER system 
ijOchsenbein et akl I2000|) . from which the photometry 
data table was easily downloaded. Fig. El shows the 
extinction-corrected color-magnitude diagram for M12 
(top left). Stars on the BHB, marked as bold points 
in Fig. El were selected by eye. M12 is relatively metal- 
poor [Fe/H] ~ —1-4, similar to our halo star samples with 
median [Fe/H]= -1.7. 

The inset in Fig. 1161 shows the luminosity function 
of BHB stars in M12. There is considerable uncer- 
tainty in the distance modulus and the metallicity of 
M12, with values ranging from (m - M) = 14.22 ± 0.11 
for [Fe/H]= -1.14 to (m - M) = 13.96 ± 0.11 for 
[Fe/H]= —1.61. We calculate absolute magnitudes using 
the distance modulus (m — M) = 14.05 ± 0.12 appropri- 
ate for [Fe/H]= -1.4 ijHargis et alJl2004|) and extinction 
E{B - V) = 0.19 ± 0.02. The median extinction to our 
halo BHB stars, by comparison, is E(B — V) = 0.03 
mag. The M12 BHB luminosity function has median 
My = 0.75, a quarter of a magnitude brighter than our 
field BHB samples. The disagreement in median My is 
significant only at the la level, however, since the Ml 2 
distance modulus, the M12 extinction correction Ay , and 
the My (BHB) zero-point are all uncertain to ±0.1 mag. 

The shape of the M12 BHB luminosity function is 
nearly identical to the shape of our field BHB luminos- 
ity functions. We use a K-S test as before and calculate 
the likelihood that the M12 and our field BHB luminos- 
ity functions are drawn from the same distribution. Be- 
cause our goal is to compare the shapes of the luminosity 
functions, we match the median My of the observed and 
M12 samples when performing the K-S test. We find 
likelihoods ranging from 24% to 44% for our field BHB 
samples, indicating that they likely share the same dis- 
t ribution of BHB lu minosities as the M12 BHB sample. 

iPiotto et al.1 l)2002f) provide on-line data for 74 addi- 
tional globular clusters that we use for further compari- 
son. The data come from Hubble Space Telescope ^43 91^ 
and F555W imaging from which lPiotto et all l|2002fl de- 
rive de-reddened B and V magnitudes. Because Hubble 
Space Telescope targeted the centers of the globular clus- 
ters, not all of the color-magnitud e diagrams are as well 
sampl ed as the [Hargis_gt alT l)2004ll M12 data. One third 
of the Piotto et al. ( 2002) globular clusters have horizon- 
tal branches that are non-existent or too sparsely popu- 
lated to provide a meaningful comparison with our BHB 
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Fig. 16. — The extinction-corrected color-magnitude diagrams for M12 (NGC 6218; IHargis et aT1l2(30l) and 15 additional globular 
clusters IPiotto et aH20021) . BHB stars are marked as bold points. The luminosity function of the BHB stars are shown in the inset panels. 



samples. Of the remaining 51 globular clusters, we select 
the 15 with E(B — V) < 0.1 to minimize uncertainties in 
extinction. These 15 globular clusters span the range of 
metallicity -2.3 < [Fe/H] < -1.2. 

Figure El shows the extinction-corrected color- 
magnitude d i agram s for the 15 globular clusters from 
IPiotto et alJ J2002). We calculate absolute magnitudes 
using the distance moduli and extinction values given by 



IPiotto et alJ |2002). Interestingly, each of the 15 globular 
clusters has a different median BHB My than our BHB 
samples. There is no correlation with globular cluster 
metallicity. The average globular cluster median BHB 
luminosity is My = 0.7±0.25. Although a 0.3 mag differ- 
ence from our BHB samples is not formally significant, we 
expected better agreement when averaging over this set 
of globular clusters. In Paper I, we noticed a similar ~0.3 
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magnitude discrepancy between our absolute magnitudes 
and the theoretical calculations for the zero ag e hori- 
zontal branch. For example, the lDemaraue et alJ (2000) 
model discussed here (Fig. [HJ is 0.1 mag brighter than 
the PSB relation. II the zero-point ol our Mv{BHB) 
relation is in error by 0.3 mag, then our BHB stars are 
14% more distant than our current estimates. Given the 
strong dependence of BHB luminosity on color, it may be 
worth re- visiting th e PSB analys i s and their zero point. 

The shape of the iPiotto et alJ l)2002|) globular cluster 
BHB luminosity functions are in good agreement with 
our field BHB luminosity functions. The globular cluster 
BHB stars are marked as bold points in Fig. El We select 
BHB stars by color/magnitude cuts, imposing the same 
limits (B — V) < 0.24 and My < 2.1 as for our samples. 
The insets in Fig. 1161 show the luminosity functions of 
globular cluster BHB stars. We perform K-S tests on the 
shapes of the luminosity functions as before, first match- 
ing the median My of the globular cluster BHB stars to 
our BHB samples. The likelihoods that the BHB stars 
are drawn from the same distribution range from 20% to 
80%, with the exception of two globular clusters. NGC 
6229 and NGC 7078 (M15) have significant extended hor- 
izontal branches and thus a much broader distribution of 
BHB luminosities than our field BHB samples. The K-S 
test yields a 10% likelihood for NGC 6229, which would 
only allow a very marginal rejection of the null hypoth- 
esis of a common parent distribution. A 1% likelihood 
is obtained for NGC 7078, indicating that its luminosity 
distribution is not consistent with our field BHB sam- 
ples. We note that the metallicities of NGC 6229 and 
NGC 7078 are [Fe/H]=-1.43 and -2.25, respectively. 

Of the sixteen globular clusters displayed in Fig. 1161 
fourteen (88%) have luminosity functions consistent in 
shape with our field BHB luminosity functions. Thus, 
under the assumption that the My-color relation is in- 
trinsic to stars on the BHB, we find that field BHB stars 
in the halo and BHB stars in globular clusters with BHBs 
appear to share a common distribution of luminosities. 
The exception to this conclusion are globular clusters 
with significant extended BHBs; we do not see large num- 
bers of extended BHB stars with (B — V)o < —0.1 in our 
samples. In the future, it would be useful to compare 
our field BHB luminosity function with dwarf sphcroidals 
and open clusters. 

5. CONCLUSIONS 

We extend the Century Survey Galactic Halo Project 
based on a new 175 deg 2 spectroscopic survey for BHB 
stars. We make use of the 2MASS and SDSS photo- 
metric catalogs, and show that the 2MASS and SDSS 
color-selection is 38% and 50% efficient, respectively, for 
BHB star s. The 2MASS sele ction for BHB stars is 65% 
complete ijBrown et alJl2004|) . but is likely to be worse 
in the magnitude range 15 < Jo < 15.5 because of large 
photometric errors scattering BHB stars out of the nar- 
row color selection range. The SDSS completeness for 
BHB stars is also magnitude dependent and appears to 
drop to 50% in the magnitude range 15 < g' < 15.5 
because of saturation problems. 

We analyze the global properties of the original Cen- 
tury Survey, the 2MASS-selected, and SDSS-selected 
BHB stars, and find them consistent with a predomi- 
nantly halo population. The median metallicity of the 



BHB stars is [Fe/H]=— 1.7. K-S tests indicates that 
the BHB samples share a common metallicity distribu- 
tion. The velocity dispersion of the BHB stars is 108 
km s _1 . The mean Galactic rotation of the BHB stars 
3 < \z\ < 15 kpc is — 4±30 km s _1 . Our samples also in- 
clude a likely run-away B7 star 6 kpc below the Galactic 
plane. 

The luminosity of a BHB star is primarily tempera- 
ture (color) dependent. The shape of the My-color re- 
lation is due to th e physics of BHB stars. We show 
that the shape of the lPreston. Shectman. fc Beersl l)1991f) 
obs ervationally-deriyed My -color relation corresponds to 
the lDem araue et alJ (1200 01) theoreti cal BHB models and 
to the iKenvonfcHartmanrJ l)1995f) bolometric correc- 
tions. We derive luminosities to our field BHB stars 
under the assumption that the My-color relation is in- 
strinsic to stars on the BHB. 

The My-color and My-metallicity relations impose se- 
lection biases on a flux-limited survey. A flux-limited sur- 
vey samples hot BHB stars with {B — V)o — —0.1 to 64% 
of the depth for BHB stars with (B-V) = +0.1. Sim- 
ilarly, flux-limited survey samples metal-rich BHB stars 
with [Fe/H]= —1 to 82% of the depth for metal-poor 
BHB stars with [Fe/H]= -3. 

We calculate the luminosity function for our field 
B HB star samples us i ng the maximum-likelihood method 
of lEfstathiou et al.l l)1988fh a technique that is non- 
parametric and unbiased by density inhomogeneities. 
The luminosity function for field BHB stars is charac- 
terized by a steep rise at bright luminosities, a peak be- 
tween 0.8 < My < 1.0, and a tail at faint luminosities. 
We show that the luminosity function is not determined 
by the shape of the My-color relation, but rather the 
way this relation is populated. We compare our lumi- 
nosity functions with the BHB luminosity functions de- 
rived from sixteen different globular clusters. K-S tests 
indicate that globular clusters with BHBs, but not globu- 
lar clusters with significant extended BHBs, have similar 
distributions of BHB star luminosities as our field BHB 
star samples. 

We plan to analyze our samples of BHB stars for ve- 
locity and spatial sub-structure. Knowing the global 
properties and luminosity function of the BHB stars is 
an important step in this analysis. Furthermore, know- 
ing the 2MASS and SDSS color-selection efficiencies and 
completenesses for BHB stars are guides our continuing 
observations. The eventual goal of our Galactic Halo 
Project is to identify star streams in the halo and thus 
to test the hierarchical picture for galaxy formation. 
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TABLE A2 
PHOTOMETRY 









Jo 


g'o 


E(B-V) 


WO 




ID 


a J2000 


<5j2000 


(mag) 


(mag) 


(mag) 


(mag) 


BHB 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



CHSS 1598 


3:41:13.2 


0:48:37 








15.27 


± 


0.013 


0.09 


-0.03 





CHSS 1599 


3:43:57.6 


0:08:57 








15.11 


± 


0.020 


0.09 


0.15 





CHSS 1600 


23:00:20.9 


-0:17:10 


14.22 


± 


0.026 








0.05 


0.16 





CHSS 1601 


23:02:10.7 


-1:01:10 


14.65 


± 


0.038 








0.05 


-0.01 


1 


CHSS 1602 


23:03:58.3 


-1:08:12 


13.59 


± 


0.030 








0.04 


-0.02 


1 



Note. — Table lA~2l is presented in its entirety in the electronic edition of the Astronomical 
Journal. A portion is shown here for guidance and content. 



TABLE A3 

SPECTROSCOPIC AND STELLAR PARAMETERS 











^radial 




T eff 


logs 




Dist 


My 


ID 


KP 


HP2 


CP 


(km s" 1 ) 


Type 


(K) 


(cm s -2 ) 


[Fe/H] 


(kpc) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 



CHSS 1598 


1.97 


11.00 


0.89 


41.0 ±27.6 


22.0 ±1.0 


8413 


4.99 


-0.17 


5.73 


1.28 


CHSS 1599 


1.72 


10.11 


1.38 


-13.9 ±26.1 


31.6 ±2.6 


8243 


4.99 


-0.53 


1.93 


3.60 


CHSS 1600 


1.90 


10.07 


0.81 


-44.4 ±11.0 


21.2 ±1.0 


8208 


4.99 


-0.29 


2.75 


2.09 


CHSS 1601 


0.47 


10.91 


0.05 


-92.4 ±10.0 


21.6 ±1.2 


9111 


3.50 


-1.77 


4.67 


1.23 


CHSS 1602 


0.31 


10.91 


0.27 


9.9 ± 9.8 


21.6 ±1.2 


9095 


3.50 


-1.49 


2.65 


1.33 



Note. — Table lA~3l is presented in its entirety in the electronic edition of the Astronomical Journal. A 
portion is shown here for guidance and content. 
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APPENDIX 
DATA TABLES 

Tables lA"2l and 1X51 list the photometric and spectroscopic measurements for the 2MASS-selected and SDSS-selected 
samples. The tables contain 257 entries and include every 2MASS- and SDSS-selected object except for the 27 G-type 
stars in the SDSS-selected sa mple . T he S DSS-selected G-types have erroneous photometry, likely due to saturation 
problems in the SDSS. Tables lA"2l and lA"3l are presented in their entirety in the electronic edition of the Astronomical 
Journal. A p ortion of the tables are shown here for guidance regarding their format and content. 

Table lA"2l summarizes the photometry. Column (1) is our identifier. The designation CHSS stands for Century 
Halo Star Survey and is chosen to be unique from previous surveys. Column (2) is the J2000 right ascension in 
hours, minutes, and seconds. Column (3) is the J2000 declination in degrees, arcminutes, and arcseconds. Col- 
umn (4) is the 2MASS extinction-corrected Jo magnitude for the 2MASS-selected stars. Column (5) is the SDSS 
extinction-corrected q' n magni t ude f or the SDSS-selected stars. Column (6) is the E(B — V) reddening value from 
ISchlegel. Finkbeiner. fc Davisl (119981). Co lumn (7) is the BVO color predicted from 2MASS or SDSS photometry and 
Balmer line strengths ( Br own et alJl200"3T) . Column (8) is the BHB classification: 1 if the star is BHB, if it is not. 

Table lA"3l summarizes the spectroscopic and stellar parameters. Column (1) is our identifier. Column (2) is the 
KP (Ca II) index. Column (3) is the HP2 (H5) index. Column (4) is the CP (G-band) index. Column (5) is the 
heliocentric radial velocity in km s _1 . Column (6) is the spectral type, where B0=10, A0=20, F0=30, and so forth. 
Column (7) is the effective temperature in K. Column (8) is the log base 10 of the surface gravity in cm s~ 2 . Column 
(9) is the metallicity given as the logarithmic [Fe/H] ratio relative to the Sun. Column (10) is the estimated distance 
in kpc. Column (11) is the absolute My magnitude corrected for reddening, given the estimated distance. 
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